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(Dated: 7th February 2008) 

The polarization of the W boson in t —> Wb decay is unambiguously predicted by the Standard 
Model of electroweak interactions and is a powerful test of our understanding of the tbW vertex. 
We measure this polarization from the invariant mass of the b quark from t Wb and the lepton 
from W ^ Iv whose momenta measure the W decay angle and direction of motion, respectively. In 
this paper we present a measurement of the decay rate (/v+a) of the W produced from the decay 
of the top quark in the hypothesis of V+A structure of the tWb vertex. We find no evidence for the 
non-standard V+A vertex and set a limit on fv+A < 0.80 at 95% confidence level. By combining 
this result with a complementary observable in the same data, we assign a limit on fv+A < 0.61 
at 95% CL. This corresponds to a constraint on the right-handed helicity component of the W 
polarization of /+ < 0.18 at 95% CL. This limit is the first significant direct constraint on fv+A in 
top decay. 

PACS numbers: 14.65.Ha, 12.15.Ji, 12.60.Cn, 13.88.-l-e FERMILAB-PUB-04-353 



The large value of the top quark mass has led to spec- 
ulation that the top quark could play a role in the mech- 
anism of the electroweak symmetry breaking Q. If so, 
the electroweak interactions of the top quark could be 
modified 0. Such a modification could alter the V— A 
structure of the tbW interaction which in turn would lead 
to an altered W polarization in top decay 0, 0, • Pos- 
sible scenarios that would introduce a V-l-A contribution 
to the tbW vertex include SU{2)l x SU{2)j^ extensions of 
the standard model . One such model invokes new mir- 
ror particles to assist a top-condensate in breaking elec- 
troweak symmetry 3. The theory of "beautiful mirror" 
fermions predicts a fourth generation up- type quark with 
right-handed weak interactions which could contaminate 
the top sample or induce a right-handed top electroweak 
interaction by mixing with the top quark 8]. 

Indirect limits of right-handed t — > bW currents have 
been placed using the process b — > 57, which proceeds via 
an electroweak radiative penguin process Q . These lim- 
its are stringent, but scenarios can be envisaged where 
other contributions to 6 — > 57 might invalidate these 
bounds. The goal of this study is a direct measurement 
of the tbW vertex from the electroweak decay of top. 

The spin-one W has three possible helicities; for the 
we label these as —1 (left-handed), (longitudinal), 
and +1 (right-handed), with the opposite convention for 
the . Because Mt > M^y, a large fraction of the 
W bosons produced in top decay will be longitudinally 
polarized The fraction is given by 



^^0 



[Mi/M^ + 2) ■ 



(1) 



For the current values of Mt = 174.3 ±5.1 GeV and 
Mw = 80.425 ±0.038 GeV [l^, this corresponds to Fq = 
0.70 ± 0.01. If there were a non-standard model V+A 
contribution to the top decay vertex, such contribution 
would not decrease the branching ratio to longitudinal W 
bosons but would instead decrease the branching ratio to 
left-handed W bosons, replacing some of this rate with 
an enhanced right-handed component. 

Leptons from the decay of longitudinally polarized W 




10000 15000 20000 



Figure 1: The theoretical distributions of Mfi, for purely V— A 
and V+A hypotheses, using the correct lepton-6 pairing. The 
M|j can be used to discriminate between the two hypotheses 
as it peaks at higher values for V+A. This ideal case does not 
include detector and trigger effects or the intrinsic lepton-6 
mass resolution. 



bosons have a symmetric angular distribution of the form 
1 — (cos-^^)^, where -0? is defined as the angle in the W 
rest frame between the lepton and the boost vector (/3) 
from the top rest frame to the W rest frame. Maxi- 
mal parity violation in the V— A electroweak theory pre- 
dicts that the non-longitudinal W helicity is purely left- 
handed in the limit of massless final state fermions. This 
creates an asymmetric angular distribution of the form 
(1 — cos'f/'^)^ [3- Due to angular momentum conserva- 
tion, even though the massive top quark may be left- 
or right-handed, positively polarized bosons are not 
possible since a massless b quark must be left-handed. 
A small right-handed component (0.04 %) of the form 
(1 + cosV'f)^ results when the mass of the b quark is 
considered. 

This analysis exploits the relationship between the an- 
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gle '0^ and the invariant mass of the £b pair, produced 
in the top decay chain t Wb, W ^ iv to determine 
the polarization of the W boson. The angle ipg can be 
related to the £b invariant mass by 

Ml = \{M^ - M2,)(1 + cosV-;). (2) 

In the V— A theory, the lepton and b jet in the W rest 
frame tend to move in the same direction, but in a V+A 
decay, the lepton and b jet typically move in opposite di- 
rections. Therefore, M^^^ would be larger on average from 
a V+A contribution as shown in Fig. ^ This difference 
can be used to determine fv+A, the fraction of t quarks 
which decay with a V+A interaction. 

If the interaction has both V— A and V+A contri- 
butions, the total angular distribution will be approxi- 
mately described by summing over weighted linear com- 
binations of the above angular distributions. The sum- 
ming of rates correctly describes the angular distribution 
from longitudinal and either a pure V+A or V— A distri- 
bution; however, if there is a combination of V— A and 
V+A interactions, they may interfere with some relative 
phase. The present analysis neglects this interference, 
which would have the largest impact for fv+A — 0.5. 
These interference effects are only of order l/7b, the 
boost of the b quark in the top rest frame, and there- 
fore are estimated to affect the angular distributions 0| 
at no more than the 10% level. The associated uncer- 
tainty is therefore not significant compared to expected 
statistical and systematic uncertainties. 

Experimentally, is a reliable observable in tt decay 
at a hadron collider because no information about the top 
or W rest frames is required, and therefore the unknown 
boost of the ti system along the beam direction does not 
disrupt the measurement. This technique also avoids the 
need to rely on the missing transverse energy (j^t) due 
to the neutrino. Thc-^T is poorly measured compared to 
other kinematic quantities in the event and is ambiguous 
in events with two final state neutrinos, e.g., both 
and W~ from the tt decay leptonically. 

The present study uses data from pp collisions at 
^fs = 1.8 TeV collected by the ColHder Detector at Fer- 
milab (CDF),1J|1 during the period 1992-1995 (Run I). 
The integrated luminosity of the data sample is 109 ± 7 
pb~^. Events were selected El HI and assigned to three 
different tt subsamples chosen for their low background 
and high efficiency for b jet identification. Each sample 
is classified by the number of leptons and identified 6 jets 
in the final state. 

The "dilepton" sample is dominated by ti in which 
both W bosons decay to an electron or muon and neutri- 
nos. Events are selected by requiring > 25 GeV, one 
muon and one electron of opposite charge with Pt > 20 
GeV in the central pseudo-rapidity region (1 7/ 1 < 1.0) 1^, 
and two jets with Et > 10 GeV and \ri\ < 2.0. This is 



a subsample of the dilepton events used in other anal- 
yses considering only e + /i + jets events in order 
to remove the dominant background, which is Drell-Yan 
production of ee or fj,fj.. The significant remaining back- 
grounds are decays to electron and muon of Z ^ rr, 
WW in association with extra jets, and W production as- 
sociated with three or more jets, where one jet is misiden- 
tified as an electron or a muon. No attempt is made to 
identify 6 jets explicitly. However, initial and final state 
gluon radiation can result in extra jets, so the b jets are 
assumed to be the two highest Et jets, which is correct 
in ~80 % of dilepton events. There are four Alei, combi- 
nations in each dilepton event. 

The other two samples used in the analysis require only 
one W to decay into an electron or muon and a neutrino 
and the other W to decay hadronically ("lepton+jets" ). 
These events are selected by requiring one electron or 
muon with Pt > 20 GeV, in the central region as above. 
At least four jets are required, three of which must have 
Et > 15 GeV, |?7| < 2.0, and the fourth must have Et > 
8 GeV and I77I < 2.4. The background for these events 
consists predominantly of direct production of a W plus 
extra jets and its behavior is modeled with the VECBOS 
generator To reduce the background, at least one 

jet must be identified as a & candidate (^-tagged) with 
a topological algorithm requiring tracks in the jet recon- 
structed with the silicon vertex (SVX) detector to form 
a secondary vertex 0, 0| . This requirement is 48 % 
efficient for tagging at least one b jet in a tt event [l8l |. 
Without any &-tag, the expected signal to background 
ratio [S/B) of the sample is 0.4, whereas requiring one 
&-tag improves S/B to 5.3. The &-tag also selects the 
jet to be paired with the lepton to form Ma,. Events 
with a single 6-tagged jet comprise the "single-tagged" 
sample, and have one measured Mu, which is correct half 
the time. Events with both b quarks tagged make up 
the "double-tagged" sample, have a S/B of 24, and pro- 
vide two Ma pairings, at least one of which combines the 
wrong b with the £. 

A total of 7 events were found in the dilepton e/i sam- 
ple with an expected background of 0.76 ± 0.21 events. 
In the single-tagged sample 15 events were found with 
a background 2.0 ± 0.7, and in the double-tagged sam- 
ple there were 5 events with a 0.2 ± 0.2 background. 
Note that since right-handed leptons have higher Pt, an 
increase in events passing the lepton Pt trigger require- 
ment could also indicate a V+A theory. However, any 
potential observed rate increase would be deemed to be a 
posteriori knowledge from the point of view of this analy- 
sis, and therefore only the shape of the A/|j distributions 
is considered. 

The M'^j^ distributions of the data are fit to a linear 
combination of three predicted Mf^^ distributions: tt pro- 
duction with a V— A interaction, ti production with a 
V+A interaction, and background. The fit maximizes a 
binned likelihood as a function of fv+A- Likelihood scans 
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Systematic Uncertainties 


1 ("\T~\ TTl fl QQ 




Tpf PTiPrP'v scalp 


0.04 


R^i p K ffTm inn '^npmp 


0.05 


rsfi p KCTTii inrl nnnn Fi M i"i nn 


0.05 


ISR gluon radiation 


0.04 


FSR gluon radiation 


0.03 


B tagging efficiency 


0.03 


Parton distribution functions 


0.02 


Monte Carlo statistics 


0.01 


Relative acceptance 


0.005 


Total systematic 


0.21 



Table I: Summary of systematic uncertainties in terms of the 
shift in measurement of the V+A fraction. The systematic 
uncertainties shown for the top mass and jet energy scale are 
after considering the correlations between the two; without 
these corrections the systematic uncertainties are 0.21 and 
0.14, respectively. 

are performed both inside and outside the physical region 
of [0, 1] in fv+A, and the level of backgrounds in each fit 
is allowed to vary within the estimated uncertainties. 

The predicted distributions are calculated sepa- 
rately for dilepton, single-tagged, and double-tagged data 
samples, by Monte Carlo simulations of tt and back- 
ground. The effects of predicted kinematics, decay dis- 
tributions, detector acceptance, and resolution are all 
considered. The HERWIG event generator [T5'| with the 
MRST h-g PDF set ^ was used to model tt production. 

For cases with two possible b jets that can be matched 
to a lepton (the dilepton and double-tagged samples), 
the fit is performed to two-dimensional distributions of 
Mff^ and Mff^ , thus taking into account that only one 
can be correct. Naively, this ambiguity in assignments 
of leptons and b quarks to one top quark would appear 
to be problematic in this measurement. However, while 
correct pairings are limited kinematically by — 
for a massless b quark, incorrect pairings often have sig- 
nificantly higher mass. With our two dimensional fit, 
mispairings only increase the statistical uncertainty in 
the fit by only 15%. 

Systematic uncertainties in the measurement enter the 
analysis primarily through the prediction of the Mn, dis- 
tributions, and are evaluated by changing assumptions 
in the Monte Carlo simulation. Listed individually in 
Table ^ all systematic uncertainties added in quadra- 
ture represent a 0.21 uncertainty in fv+A- The largest 
systematic uncertainties are from the top mass and the 
jet energy scale. Increasing the top mass will increase 
Ma in top decay. The measured uncertainty of the 
top quark mass is 5.1 GeV j^J], and an increase in top 
mass by one standard deviation increases fv+A by 0.19. 
Sources of systematic uncertainty in the jet energy scale 
include the calibration of the calorimeter, the simulation 
of the calorimeter response and the modeling of fragmen- 



tation [l^]. An increase in the overall jet energy scale 
by one standard deviation would increase fv+A by 0.14. 
However, the CDF jet energy scale has a large effect on 
the world average top mass measurement. Accounting 
for the correlation between these two effects results in a 
reduction of the systematic from jet energy scale to 0.04. 

Smaller sources of systematic uncertainties were stud- 
ied in this measurement by observing the effect in sim- 
ulated pseudoexperiments. Hard gluon bremsstrahlung 
either in the initial or final state can cause significant 
mismeasurement of the b quark jet or can produce a jet 
which can be mistaken for the b quark jet itself. The size 
of the effect was conservatively estimated by removing all 
such events from the sample in a simulated measurement. 
For samples where SVX topological 6 tagging was used, 
the effect of uncertainties in b tagging efficiency as a func- 
tion of b jet Et were evaluated. Estimated background 
rates and distributions in M|j were varied as well. The 
most important of these effects is the uncertainty in the 
mean used in the VECBOS simulation of the W^-|-iets 
background as discussed in Ref. A set of CTEQ 
and MRST ^ Parton Distribution Functions (PDFs) 
were compared to the standard PDF set of MRST h- 
g and found to cause a small spread in the measured 
fv+A- Systematic uncertainty due to the limited size of 
the Monte Carlo simulation samples is also included. 

The data and expected Standard Model distributions 
are shown for each of the three samples in Fig. El We can 
combine the statistical likelihood as a function of fv+A 
for each sample into the joint likelihood shown in Fig. |31 
The combined result for fv+A and its Icr uncertainties 
are 

fv+A = -0.2lt°:g(stat. ) ± 0.21(syst. ) (3) 

The central value depends on the true top mass, 
fv+A{Mt) = -0.21 -I- 0.037(Mt - 174.3 GeV), and the 
top mass uncertainty is reflected in the systematic error. 
This central value lies in an unphysical region, but is 
more consistent with a Standard Model V— A interaction 
for the tbW vertex than a V-l-A interaction. We can place 
a one-sided upper limit on the fraction of rate due to a 
V+A component by construction of a Neyman confidence 
band in the variable fv+A |23| . This procedure results in 
an upper limit on fv+A of 0.80 at 95% confidence level. 
With the assumption of a standard model longitudinal 
helicity fraction, this corresponds to /+ < 0.24 at 95% 
confidence level. 

W polarization in top decays has also been studied at 
CDF in the same data sample using the lepton Pt [131 ^ 
the observable to discriminate between left-handed and 
right-handed W bosons, under the assumption of a fixed 
longitudinal helicity. These two results have different se- 
lection criteria, but share largely overlapping data sets. 
In addition, the observables themselves are weakly corre- 
lated, and a large fraction of the systematic uncertainties 
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Results of SVX double tagged sample Results of SVX single b-tagged sample 
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Figure 2: Data and Standard Model Monte Carlo distributions for each sample. The last bin includes combinations greater 
than 30,000 GeV^, which are predominantly the result of incorrect pairings. Errors are statistical only. 



Combined Result: fv^^= -0.21 
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Figure 3: —2 log £ as a function of fv+A for all samples and 
for the combined likelihood fit. The result for the dilepton 
sample is fv+A ~ O.OSI'^q'^ji the single-tagged sample is 
fv+A ~ -1.91^*0^48, and for the double-tagged sample is fv+A 
= 0.63l2'ji- Errors are statistical only. 



are common. Nevertheless, the overall statistical correla- 
tion of the two results is only about 0.4. Under the sim- 
plifying assumption of Gaussian uncertainties, the com- 
bined measurement using both the Ma and lepton Pt 
approaches is that the fraction of W bosons produced in 
a V+A interaction is 

fv+A = -0.07 ± 0.37(stat. ® syst. ). (4) 

The combined upper hmit is fv+A < 0.61 at 95% con- 



fidence level. In terms of the right-handed heUcity frac- 
tion, this corresponds to /+ < 0.18 at 95% confidence 
level. The combined result is inconsistent with a pure 
V-l-A theory at a confidence level equivalent to the prob- 
ability of a 2.7(7 Gaussian statistical fluctation. 

In conclusion, we have used the measurement of Met 
in ti events to measure the polarization of W bosons in 
top decay. The results are consistent with the V— A the- 
ory of the weak interaction. The data are used to set a 
limit on the fraction of top decays mediated by a V-l-A 
interaction. This is the first result providing significant 
direct evidence against a pure V-l-A theory of weak inter- 
actions in top decay; it also provides the first significant 
limits on partial admixtures of a V-l-A interaction with 
the expected V— A reaction. 

We thank the Fermilab staff and the technical staffs 
of the participating institutions for their vital contribu- 
tions. This work was supported by the U.S. Depart- 
ment of Energy and National Science Foundation; the 
Italian Istituto Nazionale di Fisica Nucleare; the Min- 
istry of Education, Culture, Sports, Science and Tech- 
nology of Japan; the Natural Sciences and Engineering 
Research Council of Canada; the National Science Coun- 
cil of the Republic of China; the Swiss National Science 
Foundation; the A. P. Sloan Foundation; the Research 
Corporation; the Bundesministerium fuer Bildung und 
Forschung, Germany; the Korean Science and Engineer- 
ing Foundation and the Korean Research Foundation; the 
Particle Physics and Astronomy Research Council and 
the Royal Society, UK; the Russian Foundation for Ba- 
sic Research; the Comision Interministerial de Ciencia y 
Tecnologia, Spain; work supported in part by the Euro- 
pean Community's Human Potential Programme under 
contract HPRN-CT-20002, Probe for New Physics; and 
this work was supported by Research Fund of Istanbul 



7 



University Project No. 1755/21122001. 



[1] C. T. Hill, Phys. Lett. B 266, 419 (1991). 
[2] R. D. Peccei and X. Zhang, Nucl. Phys. B 337, 269 
(1990). 

[3] G. L. Kane, G. A. Ladinsky and C. P. Yuan, Phys. Rev. 

D 45, 124 (1992). 
[4] M. Jezabek and J. H. Kuhn, Phys. Lett. B 329, 317 

(1994) 

[5] C. A. Nelson, B. T. Kress, M. Lopes and T. P. McCauley, 

Phys. Rev. D 56, 5928 (1997). 
[6] For a review of V+A theories, see T. D. Lee and 

C. N. Yang, Phys. Rev. 104, 254 (1956), J. C. Rati and 

A. Salam, Phys. Rev. D 10, 275 (1974), J. Maalampi 

and M. Roos, Phys. Rept. 186, 53 (1990), R. Foot, Phys. 

Lett. B 420, 333 (1998), S. h. Nam, Phys. Rev. D 66, 

055008 (2002), Q. Shafi and Z. Tavartkiladze, Phys. Rev. 

D 66, 115002 (2002), H. S. Goh, R. N. Mohapatra and 

S. P. Ng, Phys. Lett. B 570, 215 (2003). 
[7] Specifically, G. Triantaphyllou, J. Phys. G 26, 99 (2000), 

M. Lindner and G. Triantaphyllou, Phys. Lett. B 430, 

303 (1998). 

[8] D. Choudhury, T. M. Tait and C. E. Wagner, Phys. Rev. 

D 65, 053002 (2002). 
[9] K. Fujikawa and A. Yamada, Phys. Rev. D 49, 5890 

(1994). 

[10] K. Hagiwara et al., Phys. Rev. D 66, 010001 (2002) and 
2003 off-year partial update for the 2004 edition available 
on the PDG WWW pages (URL: http://pdg.lbl.gov/ ) 

[11] Private communication with T. Tait. 



[12] F. Abe et al, Nucl. Instr. Meth. Phys. Res. A 271, 387 
(1988); D. Amidei et al, ibid. 350, 73 (1994); P. Azzi et 
al., ibid. 360, 137 (1995). 

[13] F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 80, 
2767 (1998). 

[14] F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 80, 
2779 (1998). 

[15] In the CDF coordinate system, and (f) are the polar and 
azimuthal angles, respectively, with respect to the proton 
beam direction which defines the z axis. The pseudora- 
pidity r] is defined as — In tan | . 

[16] F. A. Berends, H. Kuijf, B. Tausk and W. T. Giele, Nucl. 
Phys. B 357, 32 (1991). 

[17] F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 74, 
2626 (1995). 

[18] T. Affolder et al. [CDF Collaboration], Phys. Rev. D 63, 
032003 (2001). 

[19] G. Corcella et al, JHEP 0101, 010 (2001). Both tt sam- 
ples were generated with HERWIG, the V-f A sample us- 
ing a custom version with adjustable W helicity ampli- 
tudes. 

[20] A. D. Martin, R. G. Roberts, W. J. Stirling and 

R. S. Thorne, Eur. Phys. J. C 4, 463 (1998). 
[21] L. Demortier, R. Hall, R. Hughes, B. Klima, R. Roser and 

M. Strovink [The Top Averaging Group Collaboration], 

FERMILAB-TM-2084. 
[22] H. L. Lai et al. [CTEQ Collaboration], Eur. Phys. J. C 

12, 375 (2000). 

[23] J. Neyman, Phil. Trans. Royal Soc. London, Series A, 

236, 333 (1937). 
[24] T. Affolder et al. [CDF Collaboration], Phys. Rev. Lett. 

84, 216 (2000). 



